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ABSTRACT 

In  this  thesis,  the  algorithms  for  an  A-7E  aircraft  perfor- 
mance calculator  were  developed  and  then  implemented  on  three 
small  data  processors  of  different  programming  levels  and  storage 
capabilities . 

The  utility  of  data  is  a function  of  several  variables  in- 
cluding accuracy  and  availability.  The  problem  of  retrieving 
performance  data  from  the  Naval  Air  Training  and  Operating  Pro- 
cedures Standardization  (NATOPS)  Manuals  is  significantly  les- 
sened by  the  devices  demonstrated  in  this  investigation.  Nine 
performance  chart  groups,  yielding  data  usually  considered 
necessary  for  flight,  were  reduced  to  a series  of  analytical 
expressions.  These  analytical  expressions  were  demonstrated  to 
reproduce  NATOPS  Manual  data  to  a high  degree  of  accuracy. 

Implementation  was  demonstrated  on  a desk  computer,  a hand 
held  calculator  and  a microprocessor. 
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I . INTRODUCTION 


The  Naval  Air  Training  and  Operating  Standardization 
(NATOPS)  Manual  is  the  official  standard  of  the  United  States 
Navy  for  "...information  on  all  aircraft  systems,  performance 
data,  and  operating  procedures  required  for  safe  and  effective 
operations."  [1] 

The  purpose  of  this  thesis  was  to  develop  algorithms  of 
the  more  often  used  NATOPS  performance  charts  for  the  A-7E 
aircraft,  examine  their  accuracy  and  implement  them  on  small 
data  processors  that  might  be  adaptable  to  shipboard  or  aircraft 
onboard  use.  The  interpretation  of  NATOPS  performance  charts 
is  an  error  prone  and  time  consuming  procedure  ev^en  for  experi- 
enced users.  The  need  for  a system  to  eliminate  this  laborious 
process  has  been  fully  documented  in  a thesis  completed  in  June 
1978  by  LCDR  W.M.  Siegel  [2].  In  his  investigation,  LODR  Siegel 
devised  an  efficient  procedure  to  develop  algorithms  from  the 
NATOPS  performance  charts  and  exercised  this  procedure  on  the 
problems  of  "Takeoff  Ground  Roll  Distance"  and  "Takeoff  Airspeed". 

This  investigation  is  an  extension  of  the  aforementioned 
work.  The  original  scope  of  this  investigation  was  to  develop 
algorithms  for  eleven  of  the  most  often  used  perfoi'mance  problem 
chart  groups  and  implement  them  on  the  Texas  Instruments-59 
(TI-59)  hand  held  calculator  (HHC).  All  of  the  NATOPS  perfor- 
mance charts  were  not  reduced  because  of  research  time 
limitations.  Of  the  eleven  performance  chart  groups  studied. 
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two  performance  problems,  "Time  to  Climb"  and  "Fuel  Required 
to  Climb"  were  rejected  because  of  implementation  difficulties 
on  the  TI-59  HHC  Cdiscussed  fully  in  "Development  Difficulties"). 
Therefore,  nine  performance  chart  groups  were  reduced  to  analy- 
tical expressions  and  implemented  on  the  TI-59  HHC.  To  show 
further  possibilities  and  feasibility  of  implementation  of  the 
algorithms,  they  were  1)  fully  implemented  on  the  Hewlett 
Packard-9830  (HP-9830)  desk  computer,  2)  demonstrated  on  a 
microprocessor  (INTEL  Corporation  Microcomputer  System-48), 
and  3)  considered  for  implementation  on  the  A-7E  onboard  digi- 
tal computer  and  a microprocessor  utilizing  a recently  developed 
number  processing  chip  by  the  National  Semiconductor  Corporation 
(MM57109) . 
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II.  DEVELOPMENT 


A.  GUIDELINES 

The  scope  of  this  investigation  was  established  after  a 
firm  set  of  guidelines  was  defined. 

Being  the  official  United  States  Navy  standard  for  the 
A-7Z  aircraft,  the  .^-7E  NATOPS  Manual  was  the  sole  source  of 
performance  data  used  to  develop  the  algorithms.  As  such, 
and  being  subject  to  changes  during  the  aircraft's  life  cycle, 
the  need  for  possible  future  updates  to  the  algorithms  was 
acknowledged.  The  effective  date  of  the  NATOPS  Manual  from 
which  these  algorithms  were  developed  is  March  1975.  Since  the 
performance  data  yielded  by  the  algorithms  was  identical  to 
NATOPS  Manual  performance  curves,  the  same  restrictions  and 
limitations  apply.  For  example,  takeoff  airspeed  calculation 
restricts  the  NATOPS  Manual  user  to  trailing  edge  flap  positions 
between  20  and  40  degrees  down  (Figure  1).  For  that  reason,  one 
could  not  expect  to  calculate  the  flaps  up  takeoff  airspeed 
using  the  developed  algorithms.  An  additional  feature  provided 
by  the  algorithms  was  higher  order  interpolation.  While  the 
inexperienced  NATOPS  user  might  attempt  to  interpolate  linearly 
between  non-linearly  spaced  curves,  the  algorithms  do  not. 

An  important  guideline  for  the  user's  benefit  was  to  ensure 
the  execution  of  these  algorithms  after  ijnplementation  was 
simple  enough  so  very  little  training  was  required  for  the 
users.  Intended  users  were  Naval  Flight  Officers  and  Aviators. 
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Not  included  in  the  scope  of  this  thesis  is  an  introduc- 
tion to  the  TI-59  HHC,  HP-9830  desk  computer  and  the  INTEL 
Microcomputer  System-US;  however,  to  follow  the  computer  pro- 
grams written  for  these  devices  would  required  their  basic 
understanding. 

Another  guideline  established  was  that  the  performance 
calculators  be  light  and  small  enough  to  be  physically  suited 
for  its  environment.  For  example,  the  TI-59  calculator  and 
microprocessor  could  be  used  in  a cockpit,  briefing  room  or 
Air  Operations  Center.  The  HP-9830  desk  computer  would  be 
restricted  from  cockpit  use. 

Reliability  was  a necessary  guideline. 

To  make  algorithm  implementation  on  the  TI-59  HHC  feasible 
and  since  the  program  storing  chip,  the  Continuous  Read  Only 
Memory  (CROM),  was  limited  to  5000  calculator  program  steps, 
the  library  of  nine  programs  was  required  to  fit  into  that 
space  [3]. 

Finally,  accuracy  was  a necessary  consideration.  The 
results  obtained  from  the  algorithms  were  required  to  be  at  j 

least  as  accurate  as  following  the  performance  charts  manually. 

These  accuracy  requirements  established  were:  One  knot  of  air- 
speed, 100  feet  of  altitude  or  ground  roll  distance,  100  pounds 
of  weight,  ten  seconds  of  time  and  one  nautical  mile  of  dis- 
tance . 
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B.  PERFORMANCE  CHART  REDUCTION 


The  reduction  of  the  NATOPS  Manual  performance  curves  into 


analytical  expressions  was  accomplished  by  a historically  proven 
mathematical  procedure,  "least  squares  curve  fitting".  This 
method  was  applied  to  certain  A-7E  performance  data  by  LCDR 
W.M.  Siegel  (see  Introduction,  Section  I).  His  brief  explana- 
tion of  the  "Least  Squares  Fit  Approximation  (LSFA)"  is  included 
in  Appendix  A. 

Many  performance  charts  from  the  NATOPS  Manual  contain 
three  variables  (two  independent,  one  dependent)  and  are  de- 


picted as  a two-dimensional  space  with  the  third  dimension 


illustrated  by  a family  of  curves.  The  reduction  of  such  a 


chart  can  be  accomplished  as  follows 


1.  Determine  order  of  curves  in  family  (i.e,  second  order 


2.  Apply  LSFA  to  every  member  of  the  family  of  curves 


Since  the  order  of  the  curve  families  may  vary,  a 


general  curve  family  could  be  depicted  as  follows 


(for  curve  2,) 


(for  curve  4„) 


4.  Apply  LSFA  to  the  coefficients.  For  example,  plot  A 


respectively,  yielding 
o _r-l 


Doing  the  same  with  all  coefficients 


El! 


5. 


Given  z and  x, 


^ ^ 

^ * ®m3^*  ■*■  

y can  now  be  calculated  by; 


a.  Computing  coefficients  from  equations  generated  in 


Step  4. 

b.  Applying  coefficients  to  y = + A2X  + ...A^x^"^. 

6.  It  is  important  to  note  that  although  all  curve  family 
members  must  be  of  identical  order,  the  equations  representing 
the  coefficients  as  a function  of  "z”  need  not  be  of  similar 
order . 


Although  applying  LSFA  to  the  family  of  curves  and  then 
to  their  coefficients  was  the  normal  method  of  chart  reduction, 
it  was  not  always  used  for  the  following  reasons: 

a.  Some  charts  were  two-dimensional  (LSFA  still  used). 

b.  Some  charts  were  reduced  by  inspection. 

(1)  Linear  curve  families  with  linear  spacing. 

(2)  Time,  distance,  speed  charts  (d  = v/t). 

c.  Algorithm  anomalies  (see  "Development  Difficulties"). 

When  used,  the  LSFA  was  accomplished  by  a program  pre- 
written by  the  Hewlett  Packard  Corporation  for  use  with  the 
HP-9830.  This  program,  although  greatly  facilitating  the 
development  portion  of  this  investigation,  was  written  for  a 
two-dimensional  problem  and  had  to  be  executed  at  least  once 

for  each  curve  and  once  for  each  set  of  coefficients. 

A listing  of  all  of  the  equations  making  up  the  perfor- 
mance algorithms  are  contained  in  Appendix  C.  The  A-7E 
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performance  chart  groups  from  which  they  were  developed  are 
contained  in  Appendix  B.  They  are  in  order: 

1.  Low  Level  Cruise  Performance. 

2.  Takeoff  Ground  Roll  Distance. 

3.  Maximum  Range  Cruise  Time  and  Speed  at  Constant 
Altitude . 

4.  Maximum  Range  Cruise  Fuel  Required  at  Constant 
Altitude . 

5.  Maximum  Range  Climb  Airspeed  Schedule. 

6.  Takeoff  Airspeed. 

7.  Maximum  Refusal  Airspeed. 

8.  Optimum  Endurance  Altitude. 

9.  Cruise  Ceiling. 

Future  reference  in  this  thesis  is  made  to  algorithms  and 
programs  by  the  numbers  above. 

C.  EXAMPLE  OF  CHART  REDUCTION 

An  example  of  the  procedure  discussed  in  the  previous  sec- 
tion is  presented  below.  The  chart  chosen  for  reduction  is  the 
lower  graph  of  Figure  2,  from  Phase  II  of  the  A-7E  Cruise  Per- 
formance chart  group. 

By  inspection,  all  A^  and  A^  coefficients  are  equal  to 
zero.  The  curves  appear  parabolic  and  therefore  second  order, 
yielding  N s A^M* . The  example  follows : 

N = intermediate  result 

M s mach  number 

D = drag  count 
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Cruise  Performance 
Phase  II  NATOPS  Chart 
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CURVE  EQUATION 
N = 1.3915M* 

N s 2.7787M* 

N = 4.1658M* 

N s 5.5530M* 

N = 8.3273M* 

N = 11.102M* 

By  plotting  the  coefficients  versus  D (drag  count),  the 
LSFA  yields: 

Aj  = (4.3732E-3)  + .027743D  and  therefore, 

N = ((4.3732E-3)  + .027743D)M*. 

This  was  a particularly  simple  chart  to  reduce  but  illus- 
tlrates  the  procedure . 

D.  DEVELOPMENT  DIFFICULTIES 

The  normal  method  of  reducing  performance  curves  did  not 
always  yield  useful  information.  One  reason  was  although  the 
NATOPS  Manual  Performance  curves  were  constructed  from  experi- 
mental data,  families  of  curves  occasionally  had  very  unusual 
spacing.  They  also  were  not  always  a true  curve  family,  that 
is,  they  were  of  varying  order.  This  can  be  visually  detected 
in  the  lower  graph  of  Phase  III  of  the  A-7E  Cruise  Performance 
chart  group  (Figure  3).  The  unequal  and  varying  spacing  be- 
tween curves  with  different  "reference  numbers"  is  obvious. 
Although  the  coefficients  for  each  curve  can  be  calculated, 
the  coefficients  determined  for  a LSFA  equation  for  an  inter- 
mediate curve  would  be  incorrect.  To  be  usable  for  the  normal 
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method  of  chart  reduction,  a chart  must  have  equal,  constantly 
increasing,  or  constantly  decreasing  spacing  between  curves. 
When  such  an  incompatible  chart  was  encountered,  it  was  neces- 
sary to  interpolate  between  them.  Two  chart  groups  eliminated 
from  consideration,  "Fuel  Required"  and  "Time  to  Climb  from  Sea 
Level  to  Selected  Altitude",  contained  so  many  such  curves 
(11),  that  very  high  order  expressions  would  have  been  re- 
quired to  compute  the  coefficients,  making  implementation  on 
the  TI-59  HHC  impractical.  The  A-7E  Cruise  Performance  lower 
chart  of  Phase  I had  the  same  anomaly  (Figure  4).  Because  of 
the  importance  of  the  low  level  mission,  however,  the  algorithm 
for  this  chart  was  developed,  for  sea  level  only  though.  The 
multiple  algorithm  was  not  developed  but  could  have  been  for 
implementation  on  a desk  computer. 

Another  reason  a straight  application  of  LSFA  was  not 
always  appropriate  was  the  uniqueness  of  the  upper  graph  of 
Phase  I of  the  A-7E  Cruise  Performance  chart  group  (Figure  4). 
This  chart  requires  entry  from  the  lower  chart.  A line  is 
traced  upward  until  the  user  contacts  the  appropriate  Drag 
Count  Line  (dotted  lines).  The  first  pass  through  the  Mach 
Number  axis,  a result  of  the  lower  chart,  was  defined  M*. 
Instead  of  now  tracing  horizontally  to  the  Transfer  Scale  axis 
(this  value  defined  TS*),  one  must  trace  "between  the  solid 
guidelines"  to  the  interception  with  a line  traced  vertically 
upward  from  the  desired  Mach  number,  M.  The  Transfer  Scale 
would  now  be  manually  obtained  by  tracing  horizontally  to  the 
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vertical  axis.  To  develop  the  algorithm  for  this  problem, 
the  equations  of  the  guidelines  were  also  calculated  as  a func- 
tion of  Mach  number.  The  values  of  the  Transfer  Scale  resulting 
from  M*  intercepting  the  guidelines  and  tracing  horizontally 
to  the  vertical  axis  were  called  TS^^jTSj** 
to  bottom.  The  original  position,  (M*,  TS*),  could  now  be 
determined  in  relation  to  (M*,  and  (M*,  . "n"  and 

"n+1"  indicate  the  upper  and  lower  guidelines,  respectively, 
which  bracket  (M*,  TS*).  This  ratio  provided  the  initial  posi- 
tion relative  to  the  guidelines: 

R = (TS*-TS  ^ 'M/(TS*-TS 

n+1  n n+1 

Using  the  desired  Mach  number,  M,  the  Transfer  Scales  for 
the  same  two  enclosing  guidelines  were  calculated  (TS^  and 
TS^^^).  The  final  position  relative  to  the  guidelines  was 
maintained  using  the  original  ratio  by  solving: 

R = (x-TS  ^,)/(TS  -TS  ..)  for  x. 
n+i  n n+i 

"x"  is  the  Transfer  Scale  with  which  the  user  now  proceeds  to 
Phase  III  of  this  performance  chart  group.  Figure  5 depicts 
this  problem  graphically. 

E.  ACCURACY 

A large  number  of  results  comparisons  between  the  generated 
algorithms  and  manually  traced  performance  problems  were  made. 

An  infinite  number  of  comparisons  would  be  required  to  check 
all  possibilities,  but  since  the  mathematical  theory  was  so 
basic,  the  number  of  checks  accomplished  were  considered  suf- 
ficient . 
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All  nine  algorithms  were  checked  for  accuracy  on  the 
HP-9830  desk  computer.  The  number  of  checks  for  each  algo- 
rithm was  proportional  to  the  ease  of  manually  tracing  through 
the  performance  charts.  The  author  spent  considerable  time 

obtaining  performance  results  from  the  NATOPS  Manual  charts  ; 

and  a relatively  small  amount  of  time  computing  the  problems  | 

on  the  desk  computer  once  the  algorithms  had  been  implemented. 

In  a significant  number  of  instances,  the  results  disagreed, 
but  after  rechecking,  the  solution  obtained  manually  was  in 
error.  This  supported  the  contention  that  manual  manipulation 
of  the  performance  charts  is  an  error  prone  procedure,  even 
with  an  experienced  user. 

In  a few  rare  instances,  the  author  entered  the  required 
given  data  incorrectly  into  the  desk  computer.  These  miskeying 
errors,  not  procedural,  were  noticed  as  soon  as  the  answer  was 
produced.  A user  familiar  with  the  A-7E  performance  character- 
istics would  normally  notice  an  answer  resulting  from  grossly 
incorrect  data  input.  It  is  acknowledged,  however,  that  there 
is  no  failsafe  check  on  the  programs.  When  using  a desk  com- 
puter, the  required  input  data  can  be  printed  along  with  the 
answer  to  ensure  the  user  of  the  correctness  of  the  input  data. 

For  a hand  held  calculator,  however,  computing  a performance 
problem  twice  would  provide  a check,  which  is  what  many  NATOPS 
Manual  users  often  do.  As  with  all  computer  programs,  a desired 
result  requires  accurate  input  data. 
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Except  for  those  noted  below,  the  results  of  programs 
checked  (using  five  significant  figures)  were  indistinguish- 
able from  the  answers  obtained  by  manually  manipulating  the 
performance  charts.  Answers  produced  from  the  algorithms 
were  rounded  off  to  the  nearest  digit. 

PROGRAM  MAXIMUM  DEVIATION 

Maximum  Refusal  Speed  2 knots 

Takeoff  Airspeed  1 knot 
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III.  IMPLEMENTATION 

A.  DESK  COMPUTER 

The  use  of  a desk  computer  capable  of  producing  A-7E  per- 
formance information  within  seconds  (less  than  three  seconds 
computation  time  for  the  longest  algorithm)  would  be  ideal  for 
a squadron  briefing  room  or  Air  Operations  Center  use.  The 
HP-9830  desk  computer  was  used  for  this  implementation  stage. 

Very  little  training  would  be  required  for  personnel  to  load 
the  programs  stored  on  a cassette  tape  cartridge  and  execute 
them. 

A knowledge  of  "basic"  computer  language  is  required  to 
fully  understand  the  nine  HP-9830  programs  in  Appendix  D [4]. 

The  nine  programs  are  in  the  same  order  as  the  algorithms  of 
Appendix  C. 

Only  in  the  Low  Level  Cruise  Performance  program  are  sub- 
routines required  for  linear  interpolation  or  for  the  iterative 
method  to  find  the  Transfer  Scale  (see  "Development  Difficulties")  . 
All  other  programs  are  straight  forward,  sequential  computations. 
In  these  programs,  the  coefficients  defining  a curve  (y  = f(x)) 
for  a given  set  of  conditions  are  calculated.  That  chart 
result,  "y",  is  then  calculated  for  the  given  independent  vari- 
able "x".  The  next  chart  of  that  group  is  similarly  treated 
and  so  on  until  the  "final  result"  is  achieveo. 

The  HP-9830  programs  are  very  useful  since  they  prompt  the 
user  to  supply  the  correct  information.  Most  of  the  programs 
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"request,  then  accept"  those  inputs  required  for  the  applicable 
NATOPS  Manual  performance  chart.  The  HP-9830  then  prints  the 
data  just  entered  (ensuring  the  user  that  data  input  was  as 
desired)  followed  quickly  by  the  solution.  The  computer  is 
instantly  ready  to  receive  new  data  for  another  calculation. 

Programs  1,  2,  3,  4 and  7 (as  identified  in  "Performance 
Chart  Reduction,  section  II-B),  are  written  in  this  "request, 
then  accept"  format.  The  shorter  programs,  5,  6,  8 and  9, 
were  written  with  an  initial  set  of  input  data  already  in  the 
program.  This  format  allowed  the  computer  to  step  incrementally 
through  the  allowable  range  of  values  for  the  input  data,  thus 
calculating  a "table  of  performance  data"  for  the  applicable 
performance  chart  group.  These  programs  are  easily  altered  to 
the  "request,  then  accept"  format  by  some  simple  edit  commands 


The  variables  used  in  the  programs  are  defined  following 
each  program  in  Appendix  D. 


B.  HAND  HELD  CALCULATOR 

The  many  favorable  features  of  the  hand  held  calculator 
encouraged  its  implementation  of  the  performance  algorithms. 

Its  small  size  allowed  consideration  for  use  in  the  cockpit. 

Its  simplicity  and  reliability  was  an  advantage  making  it 
especially  suited  for  users  of  varying  experience  (including 
no  experience).  Although  its  execution  speed  was  the  slowest 
of  all  devices  used,  the  computation  time  was  still  much  faster 
than  using  the  NATOPS  Manual. 


The  Texas  Instruments-59  (TI-59)  programmable  hand  held 
calculator  (HHC)  was  selected  for  implementation.  This  selec- 
tion was  made  for  several  reasons.  At  the  time,  it  was  the 
only  calculator  available  to  the  author  which  allowed  permanent 
program  storage  (on  magnetic  cards).  Additionally,  the  Texas 
Instruments  Corporation  had  the  capability  to  combine  all  pre- 
written performance  programs,  up  to  a 5000  program  step  limit, 
onto  a Continuous  Read  Only  Memory  (CROM)  chip,  making  the 
A-7E  performance  programs  a permanent  part  of  the  calculator. 
This  CROM  chip  can  also  be  used  on  the  less  expensive  TI-58  HHC. 
These  features  made  the  TI-S8/S9  (with  CROM)  a practical  system 
for  the  .A-7E  Naval  Aviation  community. 

One  might  consider  the  calculator's  inability  to  prompt  the 
user  for  inputs  a shortcoming  of  this  implementation  candidate, 
but  a company  spokesman,  Mr.  Richard  Cuthbert , stated  a new 
face  could  be  fitted  onto  the  calculator,  identifying  different 
buttons  with  the  input  data  categories  such  as  GW  for  gross 
weight,  FLPS  for  flap  position,  T(°C)  for  temperature,  and  so 
on  [ 3 ] . 

Some  time  was  required  for  the  author  familiarization  with 
the  TI-59  HHC  and  its  capabilities.  For  a detailed  explanation 
of  comments  in  this  section  involving  TI-59  programming  and 
Appendix  E,  consult  the  user's  manual  [5]. 

All  programs  were  entered  with  the  calculator  memory  parti- 
tioned to  allow  879  program  steps  and  ten  memory  storage 
locations.  The  loss  of  program  steps  in  order  to  provide 
coefficient  storage  locations  (ten  to  one)  was  the  reason  for 
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partitioning  in  this  manner.  Only  five  significant  figures 
were  considered  necessary  for  computational  accuracy.  Consider- 
ing the  number  possibilities  (1.2345  to  1.2345E-12)  might  take 
from  six  to  ten  program  steps,  this  was  less  than  the  absolute 
ten  program  steps  sacrificed  for  a storage  location.  The  ten 
memory  storage  locations  were  used  to  store  the  input  data  at 
program  execution  start  but  were  often  reused  after  the  input 
data  storage  was  no  longer  required. 

The  programming  language  level  of  the  TI-59  HHC  is  below 
the  HP-9830 ’s  and  above  a microprocessor's  (discussed  later) 
in  sophistication.  The  algorithms  were  computed  in  a more 
space-saving  manner  than  on  the  HP-9830.  For  example,  in  com- 
puting a firsc  order  polynomial,  the  HP-9830  program  functioned 
as  follows: 

B(0)  = 

8(1)  = ^22^ 
y = 8(0)  + 8(l)x. 

The  TI-59  HHC  was  programmed  to  compute  as  follows: 

(All  * ^12^^  * ^^21 

In  the  Low  Level  Cruise  Performance  program,  the  linear 
interpolation  and  iterative  methods  to  follow  guidelines  (dis- 
cussed in  previous  section)  was  still  accomplished  using  the 
more  tedious  TI-59  HHC  language. 

Using  the  partitioning  already  described,  a program  limit 
of  879  program  steps  was  imposed  (filling  two  magnetic  cards). 
Two  programs,  "Takeoff  Ground  Roll  Distance"  and  "Low  Level 
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Cruise  Performance",  exceeded  this  limit  and  had  to  be  continued 
on  extra  cards.  These  programs  were  written  to  allow  storage 
of  an  intermediate  result  into  the  T-register.  The  rest  of  the 
cards  could  then  be  read  in,  any  lost  or  newly  acquired  input 
data  entered,  and  program  execution  would  continue,  automatically 
retrieving  the  stored  intermediate  result  from  the  T-register. 
These  artificial  necessities  for  program  completion  using  the 
magnetic  cards  would  not  be  necessary  if  the  programs  were 
stored  permanently  in  the  CROM. 

The  total  number  of  steps  required  for  the  nine  performance 
algorithms  programmed  on  the  TI-59  HHC  was  S461  steps.  By  sub- 
routining (340  steps  of  programming  are  common  to  two  programs), 
the  total  number  could  be  reduced  to  5121  steps.  The  elimina- 
tion of  the  artificial  steps  required  for  the  oversized  programs 
would  reduce  the  overage  more.  The  sole  intent  of  this  imple- 
mentation phase  was  not  to  fit  these  nine  programs  into  the 
5000  step  CROM.  If  the  inclusion  of  all  nine  programs  was 
desired,  streamlining  aid  offered  by  engineers  from  the  Texas 
Instruments  Corporation  plus  the  reduction  of  significant 
figures  in  a non-critical  area  would  accomplish  this. 

The  program  listings,  storage  location  usage,  user  instruc- 
tions, and  execution  times  are  included  in  Appendix  E. 

C.  MICROPROCESSOR 

1 . Single  Board  Computer  using  Software  for  Mathematical 
Operations  ~ ~ 

The  single  board  computer  (SBC)  implementation  was 
investigated  both  as  an  extension  of  thesis  work  and  to  meet 
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the  course  objectives  of  AE--4900  , Air  Data  Systems.  Work 
toward  this  effort  was  also  done  by  LCDR  W.M.  Siegel.  The 
performance  algorithms  were  to  be  processed  on  a SBC  using  an 
INTEL  Corporation  8048  Programmable  Read  Only  Memory  (PROM), 
external  random  access  memory  (RAM)  and  a program  counter. 
Software  development  was  completed  on  the  INTEL  Prompt-48 
(Microcomputer  System-48  language)  using  an  INTEL  8035  arith- 
metic logic  unit  (ALU).  Although  a SBC  using  the  8048  PROM 
and  requiring  a digital  keyboard  and  display  was  never  actually 
constructed  because  of  the  time  limitations,  the  software  opera- 
tion was  successfully  demonstrated  on  the  Prompt-48. 

To  preserve  the  programs  between  operation  periods , 
the  Prompt  48  was  hand  wired  as  specified  in  the  user's  manual 
to  an  ASR-35  Teletype  set  which  allowed  paper  tape  storage  [6]. 
The  Prompt-48  provided  1024  by  two  bytes  of  RAM  and  S4  by  two 
bytes  of  resident  memory.  Although  the  MCS-43  instruction  set 
will  not  be  discussed  in  this  thesis,  a basic  understanding  of 
assembly  level  language  is  necessary  to  understand  the  developed 
software  presented  in  Appendix  F [7].  This  microprocessor 
progrdun  listing  includes  the  MCS-48  instructions  in  hex  code 
and  literal  mneumonics  and  includes  full  documentation  to  facili- 
tate interpretation. 

A full  performance  algorithm  was  not  implemented  on  the 
Prompt-48  because  of  its  memory  storage  limitations.  The 
original  intent  was  to  exercise  the  software  of  the  complete 
A-7E  Takeoff  Ground  Roll  Distance  algorithm  on  the  Prompt-48. 
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After  the  necessary  routines  were  written  and  stored,  only  room 
for  three  coefficients  remained  (98  coefficients  required  for 
this  algorithm).  Since  implementation  capability  was  the 
desired  result,  the  computation  of  a second  order  polynomial 
was  considered  sufficient.  Although  this  effort  was  software 
oriented,  the  necessary  RAM  storage  for  the  additional  coeffici- 
ents and  executive  routine  could  have  been  easily  provided  for 
a SBC. 

The  software  development  for  algorithm  implementation 
required  routines  for  input/output  (I/O),  executive  direction, 
binary  to  binary  coded  decimal  (BCD)  and  BCD  to  binary  conver- 
sions, and  floating  point  binary  addition  and  multiplication 
routines.  The  I/O  and  executive  routines  were  wt'itten  by  LCDR 
Siegel.  The  nonavailability  of  a number  oriented  microprocessor 
at  the  time  of  this  effort  required  the  development  of  the 
mathematical  package  described  above.  The  advantages  for  such 
a capability  will  be  discussed  in  the  following  section. 

In  addition  to  the  microprocessor  software  developed  by 
the  author  and  LCDR  Siegel,  the  I/O  and  display  routines  would 
require  alteration  for  SBC  implementation  since  a digital  dis- 
play and  keyboard  would  replace  the  Prompt-48. 

Figure  6 illustrates  the  solution  method.  Figure  7 is 
a flow  chart  of  the  program  execution  sequence.  Figures  8 and  9 
show  the  Prompt-48  RAM  and  resident  register  memory,  respec- 
tively. 
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RANDOM  ACCESS  MEMORY  M^P 


ADDRESS 

000-069 

06A-06F 

070-079 

07A-0C6 

0C8-0C2 

0E5-0FF 

100-2EC 

300-3FF 


INPUT  AND  DISPUY 
EXECUTIVE  ROUTINE  SEGMENT 
COEFFICIENT  STORAGE 
MAIN  EXECUTIVE  ROUTINE 
BINARY  TO  BCD  EXECUTIVE  ROUTINE 
MISCELLANEOUS  SUBROUTINES 
ADDITION  AND  MULTIPLICATION  SUBROUTINES 

BCD  TO  BINARY  EXECUTIVE  ROUTINE  AND  CONVERSION  SUBROUTINES 


Figure  3 

Random  Access  Memory  Map 
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USE 
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LSB 

30 

21 

X- LOCATION 

31 

22 

ARITHMETIC  REGISTER 

32 

LSB 

23 

MSB 

33 

DISPLAY  HEX 

24 

EXPONENT 

34 

BUFFER 

25 

LSB  Y- LOCATION 

35 

^GB  1 

26 

MSB  ARITHMETIC  REGISTER 

36 

DECIMAL  POINT  MASK  | 

1 

27 

EXPONENT 

37 

CHARACTER  COUNTER  j 

1 

28 

LSB  BCD-BINARY 

38 

LSB 

29 

MSB  CONVERSION 

39 

2A 

EXPONENT 

3A 

DISPLAY  1 

1 

2B 

3B 

i 

2C 

3C 

BIT 

2D 

3D 

2E 

3E 

PATTERNS 

2F 

3F 

MSB 

Figure  9 

Resident  Register  Map 
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The  second  order  polynomial,  y = B(0)  + B(l)x  + B(2)x*, 
was  calculated  using  a mathematical  executive  routine  (alter- 
able for  any  size  polynomial  and  any  number  of  polynomials). 

The  only  mathematical  operations  required  were  multiplication 
and  addition  of  positive  or  negative  numbers.  For  speed,  binary 
arithmetic  was  used.  For  increased  storage  capability  and 
mathematical  efficiency,  a floating  point  capability  was  in- 
cluded . 

The  calculation  routine  proceeded  as  follows: 

B2*x  = (B,x) 

(3,x)  + (B, ) = (B,x  + B,  ) 

4.  i.  ^ 1 

(B^x  + 32^)*x  = (B^x^  + Bj^x) 

(B,x*  + Bj^x)  + Bq  = (B^x^  + B^^x  + Bq) 

Although  all  mathematical  operations  are  performed  in 
the  3-bit  (2-byte)  accumulator  register  of  the  30  35  .ALU  (for 
a SBC,  the  8043  PROM),  a working  accumulator  using  five  regi- 
sters (resident  memory  registers  two  through  six),  was 
established.  .All  numbers  in  the  program  (independent  variable 
"x"  after  conversion  to  binary,  coefficients  stored  in  RAM 
070-079  and  the  'result')  were  in  one  of  two  binary  conventions. 
While  in  storage,  the  numbers  were  in  "storage"  convention. 

The  numbers  were  shifted  from  "storage"  to  "working"  convention 
only  when  transferred  from  the  X and  Y locations  (see  resident 
register  memory  map.  Figure  9)  to  the  working  accumulator 
(registers  two  through  six).  When  the  desired  operation  was 
completed,  the  result  was  returned  to  the  "storage"  convention 
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and  moved  to  the  "X"  location.  Figure  10  displays  the  "storage 
and  working"  conventions . 

This  software  was  successfully  demonstrated  on  the 
Prompt-48.  The  user  instructions  for  the  Prompt-48  to  repeat 
the  demonstration  are  listed  below: 

(1)  Ensure  the  8035  ALU  or  8048  PROM  is  inserted  in  the 
"execution"  socket  of  the  Prompt-48. 

(2)  Enter  the  program  in  hex  code  in  the  proper  storage 
locations  as  listed  in  Appendix  F. 

On  the  Prompt-48,  press  the  following  keys  to  clear 
the  resident  register  memory: 

"C" 

"Registers" 

"0" 

II  It 
» 

11411 

"8" 

Do  not  press  "Program  Memory"  instead  of  "registers" 
or  the  program  just  entered  will  be  erased. 

(3)  To  execute  the  program,  press  the  following  keys: 

It  A It 

ff  2 

"Execute" 

"Go" 

"No  Break" 

"0" 

"Execute" 


I . '-Ml. Ml.  I nu  u i .«■»  — , 
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( 

f 

i 

I 

BINARY  CONVENTIONS 


Each  large  block  depicts  1 byte  which  Includes  4 bits.  The  compart- 
mented  blocks  represent  1 bit. 


STORAGE  CCTJVENTION 


In  storage  convention,  the  mantissa  is  left  justified  to  bit  5.  A 
positive  number  is  denoted  by  0 in  the  first  bit  of  the  second  byte (sign 
bit);  a 1 indicates  a negative  number. 


WORKING  CONVENTION 


In  working  convention,  the  mantissa  is  left  justified  to  bit  4. 

The  sign  bit  is  stored  in  FO(X- location  number)  and  FI (Y-location  num- 
ber) flags  of  the  program  status  word. 


The  display  will  blank,  awaiting  the  input  of  the 
independent  variable  "x".  To  enter  "x",  enter  the  digit  keys 
for  numbers  (base  10)  and  "D"  for  decimal  point.  "x"  will  be 
displayed  on  the  digital  display  as  it  is  entered.  To  compute 
the  algorithm  (second  order  polynomial),  press  "E”.  The  answer 
will  rapidly  appear.  To  calculate  the  polynomial  with  a new 
value  for  "x" , start  at  Step  3. 

(4)  To  prevent  the  time  consuming  reloading  of  the  pro- 
gram, it  is  advisable  to  store  the  program  on  a peripheral 
device  (paper  tape,  disc,  etc.). 

2 . Single  Board  Computer  using  Number  Oriented  Micropro- 
cessor 

Very  recently,  the  National  Semiconductor  Corporation 
began  production  of  a chip  intended  for  use  in  number  processing 
applications  [8].  This  chip,  the  MM57109  MOS/LSI,  is  capable 
of  all  scientific  calculator  functions,  test  and  branch  capa- 
bilities, internal  number  storage,  and  T/0  instructions.  Of 
the  specific  calculator  functions,  only  addition,  subtraction 
and  multiplication  would  be  used. 

A SBC  using  this  chip  would  need  the  8048  PROM  for 
coefficient  and  executive  routine  storage  but  would  not  need 
the  space  consuming  mathematical  package  of  the  SBC  in  the  last 
section.  A program  counter  would  still  be  required  but  external 
RAM  would  not.  The  computation  time  would  be  increased  over  the 
demonstrated  SBC  (approximate  computation  time  of  a HHC),  but 
the  simplicity  of  programming  would  make  this  proposed  SBC  very 
attractive . 
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D.  A-7E  TACTICAL  COMPUTER 


In  February  1978  the  author  made  a trip  to  the  Naval  Air 
Facility  at  China  Lake,  California.  The  purpose  of  this  visit 
was  to  receive  indoctrination  on  the  TC-2/2A  tactical  computers 
and  obtain  a programming  manual  for  these  devices.  The  desired 
goal  was  implementation  of  selected  performance  algorithms  on 
the  laboratory  bench  computer  run  by  the  A-7  Program  Office  of 
the  Naval  Weapons  Center  (NWC).  A thorough  understanding  of 
the  computer's  capabilities  and  limitations  was  provided  by  Mr. 
Robert  Westbrook,  a software  technician. 

The  A-7E  computer  provides  very  accurate  navigation  and 
weapons  guidance  capability.  The  TC-2  and  TC-2A  computers  are 
a generation  apart,  the  TC-2A  being  over  two  times  faster  and 
having  twice  the  storage  capability  of  its  earlier  version. 

Both  computers  are  operational  at  this  time.  Specific  design 
and  programming  information  is  available  from  the  programming 
manual  [9]. 

The  instruction  set  of  the  tactical  computer  provides  fixed 
point  arithmetic,  logical  transfer  of  control  (branching), 
address  modification  and  single  word  input/output  instructions 
specifically  intended  for  operations  primarily  involving  arith- 
metic. These  features  made  the  implementation  of  algorithms  a 
logical  decision.  Several  factors  made  this  implementation  by 
the  author  impractical.  The  computer  design  was  quite  old,  the 
instruction  set  being  very  tedious  and  difficult  to  interpret. 
The  computer's  inability  to  function  using  floating  point 


arithmetic  would  require  a significant  software  effort  in  that 
area  alone.  The  time  required  to  become  fully  familiar  with 
the  instruction  set,  write  the  software,  and  load  and  test  the 
programs  at  NWC  would  have  been  prohibitive  for  this  investiga- 
tion. 

It  is  hoped  that  the  programmers  at  NWC  will  be  able  to 
implement  those  algorithms  deemed  desireable  to  achieve  an 
onboard  capability.  Takeoff  Airspeed  and  Maximum  Refusal  Air- 
speed are  considered  ideal  for  implementation. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


Nine  of  the  A-7E  NATOPS  Manual  performance  chart  groups 
were  reduced  to  a series  of  analytical  expressions  or  algo- 
rithms. These  algorithms,  accurate  to  five  significant  figures, 
are  as  accurate  as  results  obtained  by  manual  manipulation  of 
the  performance  charts. 

Implementation  was  made  on  three  data  processors  of  different 
programming  levels  and  storage  capabilities.  These  devices  and 
degrees  of  implementation  were: 

(1)  HP-9a30  Desk  Computer  - complete  implementation 
with  successful  demonstration. 

(2)  TI-59  Hand  Held  Calculator  - complete  implementation 
with  successful  demonstration. 

(3)  Microprocessor  - partial  implementation  with  suc- 
cessful demonstration. 

In  view  of  the  success  of  this  investigation,  recommendations 
concerning  implementation  possibilities  are  listed  below: 

(1)  Complete  reduction  of  the  NATOPS  Manual  performance 
charts  could  be  accomplished  and  implemented  onto  a desk  com- 
puter as  one  large  program  capable  of  performance  data  computation 
within  seconds.  The  desk  computer  would  be  ideal  for  mission 
planning  on  a squadron  or  air  wing  level  or  for  Air  Operations 
Center  use. 

(2)  The  programs  written  for  the  TI-59  HHC  could  be  con- 
solidated onto  a CROM  and  used  with  a TI-58  HHC  for  use  on  a 
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squadron  level.  As  an  alternative,  the  software  could  be 
rewritten  for  any  HHC  of  comparable  capability. 

(3)  Although  implementation  on  a single  board  computer 
using  a number  oriented  microprocessor  is  completely  feasible, 
because  of  programming  ease  and  cost  consideration,  the  HHC  is 
considered  a superior  implementation  possibility  at  this  time. 

(4)  The  A-7E  tactical  computer  could  easily  be  pro- 
grammed by  software  engineers  at  NWC,  China  Lake,  California, 
to  produce  an  onboard  capability. 
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APPENDIX  A 


Least  Squares  Fit  Approximation 

References  10  and  11  describe  the  Least  Squares  Fit  Approxi- 
mation in  detail.  In  general  the  problem  is  to  represent  a set 
of  "n”  data  ’^oints  in  two-dimensional  space 

Xi*  Yi  i = 1 to  n 

by  a polynomial  expression  of  a curve  whose  degree  is  less  than 
"n".  Two  classes  of  problems  exist: 

(1)  Linearly  independent  - those  in  which  the  degree 
(d)  of  the  polynomial  is  one  less  than  the  number  of  data 
points 

d = n-1  (1) 

(2)  Linearly  dependent  - those  in  which  the  degree  (d) 
is  less  than  n-1 

d < n-1  (2) 

As  an  example,  a set  of  four  (4)  data  points  randomly  spaced 
was  chosen.  If  a third  degree  polynomial  of  the  form 

Y = A + BX  + CX*  + DX’  (3) 

were  desired,  and  the  data  points  X^  and  Y^  were  inserted 
(i  = 1 to  4)  into  four  such  equations,  an  exact  solution  for  the 
four  unknown  coefficients  would  exist.  These  four  unknowns 

could  be  found  from  the  four  equations  by  numerous  conventional  j 

techniques  (Direct  substitution,  Cramer's  rule,  etc.).  The  poly-  j 

nomial  expression  generated  would  be  termed  a "col-location"  | 

polynomial  because  its  plot  would  pass  through  all  data  points.  j 
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It  is  often  advantageous  to  describe  a set  of  data  points 
by  a curve  that  does  not  pass  through  each  point . This  type 
of  polynomial  would  be  termed  a "regression"  equation.  For  any 
set  of  data  points  an  infinite  number  of  regression  expressions 
exist  for  any  specified  degree  (except  the  linearly  independent 
case)  and  the  object  of  the  Least  Squares  Method  is  to  find  the 
polynomial  coefficients  of  the  chosen  degree  that  best  describe 
the  data  points.  In  the  previous  example  of  four  data  points, 
assume  that,  instead  of  the  third  degree  form  chosen,  a second 
degree  equation  were  selected  of  the  form 


With  four  data  points,  the  polynomial  is  overspecified  and  thus 


linearly  dependent.  For  this  case  an  infinite  number  of  solu 


a total  squared  error  term  (E)  could  then  be  defined  by  squaring 
and  summing  the  terms  attained: 


If  E were  them  minimized  for  any  given  degree  chosen,  the  best 
Least  Squares  Fit  would  have  been  achieved. 

If  the  values  for  6 from  Equation  5 were  inserted  in  Equation 
6 and  the  partial  derivative  of  E were  taken  with  respect  to  the 
coefficient  A,  an  equation  would  be  generated  that  when  set  equal 
to  zero  (0)  would  define  a minimum  value  of  E for  a given  value 
of  A.  If  the  same  operation  were  performed  with  respect  to  the 


I 


coefficients  B and  C then  three  equations  would  be  generated 
with  three  unknowns  (A,  B and  C).  The  solution  of  these  simul- 
taneous equations  would  produce  the  coefficients  A,  B and  C, 
that  would  minimize  the  value  of  E and  hence  would  produce  a 
Least  Squares  Fit  approximation  to  a set  of  linearly  dependent 
equations . 

A numerical  procedure  has  been  developed  to  accomplish  this 
task.  An  example  of  this  procedure  has  been  included  in  the 
following  paragraphs  ClO,  11]. 


Least  Squares  Fit  Method  Example 
Given  the  following  set  of  data: 

X 0 1 2 4 7 

f(X)  = Y 0 1 3 12  20 

fit  a curve  of  the  form 

f(X)  = Y - A + BX  + CX* 

STEP  1:  Substitute  all  pairs  of  data  into  the  form  equation 
yielding  the  fact  that  the  coefficients  (A,  B and  C)  must 
satisfy  all  the  following: 

0 = A + B(0)  + C(0)* 

1 s A + B(l)  + C(l)* 

3 = A + B(2)  + C(2)* 

12  = A + B(4)  + C(4)* 

20  = A + B(7)  + C(7)* 

Now  multiply  each  expression  by  its  coefficient  of  A in  that 
expression  and  add  all  equation  yielding 


ma\ 
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36  = 5A  + 14B  + 70C 

Now  multiply  each  expression  by  its  coefficient  of  B in  that 
expression  and  add  all  the  equations  yielding 

0 = 0(A)  + 0(B)  + 0(C) 

1 = A + IB  + 1C 


6 = 

2A 

+ 4(B) 

+ 8(0 

It 

CO 

4A 

+16(B) 

+ 64(0 

140  = 

7A 

+44(B) 

+ 343(0 

195  = 

14A 

+70(B) 

+ 416(0 

Now  multiply  each  expression  by  its  coefficient  C in  that 
expression  and  add  all  the  expressions  yielding 

0 = 0(A)  + 0(B)  + 0(C) 

1 = 1(A)  + KB)  + 1(C) 

12  = 4(A)  + 8(B)  +16(0 

192  =16(A)  +64(B)  +256(0 
980  =49(A)  +343(B)  + 2401(O 
1185  = 70A  + 416B  + 2674C 

Now  solve  the  following  three  previously  generated  equations  for 
the  coefficients  A,  B and  C yielding 
36  = 5A  + 14B  + 70C 
195  =14A  + 70B  +416C 
1185  =70A  +416B  +2674C 

A = -.99,  B = 2.6,  C = .065 

and 

Y = -.99  + 2.6X  + .065X* 
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The  following  plot  and  chart  depict  the  original  data  and 

I 

the  data  obtained  from  the  equation  for  the  fitted  curve:  ' 


Original 


Fitted  Curve 
Polynomial 


-.98 


1.67 


4.48 


10.46 


20.41 


Q.E.D. 


APPENDIX  B 


NATOPS  Manual  Performance  Charts 

These  charts  from  which  the  performance  algorithms  were 
developed  are  listed  below  in  order: 

Figure  Title 

B1  Cruise  Performance,  Phase  I 

B2  Cruise  Performance,  Phase  II 

B3  Cruise  Performance,  Phase  III 

B4  Cruise  Performance,  Phase  IV 

BS  Takeoff  Factor 

B6  Takeoff  Ground  Roll  Distance 

B7  Adjusted  Takeoff  Ground  Roll  Distance 

B8  Maximum  Range  Cruise  at  Constant  Altitude  (Time,  Speed) 

B9  Maximum  Range  Cruise  at  Constant  Altitude  (Fuel 

Required) 

BIO  Military  Power  Climb  Schedule 

Bll  Takeoff  Speed 

B12  Maximum  Refusal  Speed 

B13  Cruise  Ceiling  and  Optimum  Endurance  Altitude 
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NAVAm  01  4SAAE-1 


PHASE  111  - POUNDS  OF  FUEL  PER  NAUTICAL  MILE 


A-7C 

DATA  BASIS  FLIGHT  TEST 
DATE:  NOVEMBER  1971 


ENGINE;  TE41.A-2 
fUEE  E'lRAOE;  JP-S 

Fuel  density  c.s  ib/qal 


NAVAIR  OI-ASAAf-l 


TAKEOFF  FACTOR  (A-7E) 


MOOiL:  Alt 

data  AAStS  FLIGHT  TItT 
OATf  NOVfMtIA  \91\ 


INQiNI  TF4t  A-2 
PUlLCnADC  JPi 
FUCt  OCNSITV  6 • Lt/OAL. 


Il^lt  Chant*  ^ 


Figure  B5 

Takeoff  Factor 


< 
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NAVAIR  01-4SAAE-1 


TAKEOFF  GROUND  ROLL  DISTANCE  (A-7E) 


CONDITIONS; 

LSVtL  HARO  SURRACf  RUNWAV 
MILITARY  RATIO  THRUST 
LANOINO  CONIIOURATION 
SIRO  HiADWINO 
COi  aSH  MAC 
RULL  PLARt 


tNOINC;  TP4l«A*a 
PUfLGRAOl:  JP«6 
PUiL  OCNSITY  M LI/OAU 


MOOCL;  A-71 

DATA  CASIS  PLIGHT  TIST 
OATS:  NOVIMSIR  ISTI 


P«r  KumMUty  Hftctt  on  iMtoff  ditttnc*.  fround  roll 
dWnc—  inoold  mcfOMd  l\  lor  ••eh  t0%  •ncioow 
In  iht  lofotwo  liumirfitv  abcw*  40%. 


Chini*  • 


Takeoff  Ground  Roll  Distance 


NAVAIR  01-4SAAE-1 


TAKEOFF  GROUND  ROLL  DISTANCE  (A-7E) 

ADJUSTED  GROUND  ROLL  DISTANCE 

MOOCL.  A-7i  CONDITIONS:  tNGINC;  TSSI-A-} 

DATA  SASIS.  SLIONT  TtST  HARO  SUMACS  RUNWAY  SUSL  ORAOS:  r-t 

OATS:  NOVSMSSR  )97t  MILITARY  RAT60  THRUST  SUSL  OSNSITY:  S.S  LS/OAL. 

LANOINO  CONSIOURATION 
LIAOINO  SOOS  SLAPS  DOWN 


NOTI» 


Por  hymidity  «ff«ca  on  tohooff  dittmco.  ground  roll 
dIoioneM  iftould  bo  incrtawd  1%  for  MCf)  101^  inortoio 
tn  thg  rolotMO  KwffMdltv  Oout  40%. 


QROUNO  ROLL  OISTANCl  1.000  OT 


tQ94f\04  0001 
IHOliM  SSOHO  10VWIAV 


Maximum  Range  Cruise  at  Constant  Altitude  (Time,  Speed) 


FUEL  REQUIRED 


(NOINC:  TP41-A.2 
rUCL  CIIA06;  JP-5 
rU€L  OCMSITY:  6^  LB/OAL 


Maximum  Range  Cruise  at  Constant  Altitude  (Fuel  Required) 


n-4« 


Figure  BIO 

Military  Power  Climb  Schedule 
60 


L 


NAVAIR  0V4!;aAE  t 


I 


21 


Maximum  Refusal  Speed 


Figure  B13 

Cruise  Ceiling  and  Optimum  Endurance  Altitude 


APPENDIX  C 
Generated  Algorithms 

LOW  LEVEL  CRUISE  PROGRAM 

Phase  I 

Ml  = -92.512  + 236.896G 
Transfer  Scale  Versus  Drag  Count 

AO  = -2.3287  -.26316D  + .0073327D*  -(  7 . 513E-5)D^  + ( 3 . 5 396E-7 ) D'* 
-(7.78E-10)D®  + (6 .462E-13)D® 

A1  = 4.835  + 1.0956D  -.030653D*  + ( 3 . 1912E-4)D*  -( 1 . 5276E-6  )D‘^ 
+(3.408E-9)D®  -(2.8692E-12)D® 

A2  = 10.284  -1.0719D  + .031094D^  -(  3 . 2878E-4  )D^  + ( 1 . 59  5E-6  ) D** 
-(3.6009E-9)D®  + ( 3 . 0634E-12)D® 

SI  = AO  + (A1)(M1)  + (A2)(M1)^ 

Transfer  Scale  Versus  Guidelines 

BO  = 22.819  -31.7341  + 41.331^  -5.09531® 

B1  =-154.98  + 217.511  -261.731^  + 35.9051® 

B2  = 405.08  -525.561  + 607.491®  -88.7371® 

B3  = -445.62  + 542.981  -611.551®  + 92.8941® 

B4  = 184.78  -204.421  + 225.891®  -35.1891® 

S = BO  + (B1)(M1)  + (B2)(M1)®  + (B3)(M1)®  + (B4)(M1)- 

Phase  II 

R = S + 2C(4.3732E-3)  + .027743D]M® 


Phase  III 


BO  = 5.6253  -1.989R  + 3.0252R^  -1.0761R^  + .17675R‘’  -.013095R® 
+ (3.526E-4)R® 

B1  = 205.3012  -248.9317R  + 91.66355R*  -15.55218R’  + 1.224432R‘‘ 
-.0395333R®  + ( 2 . 89 6 385E-4 ) R® 

B2  = -1052.123  + 1231. 24R  -487.4233R*  + 91.6522R^  -8.662962R'* 

+ .3953974R®  - . 006905535R® 

B3  = 1680.142  -1950. 139R  + 788.8513R^  -152.5733R^  + 15.03819R'* 
-.7274139R*  + .013707R® 


R3 

= R 

R1 

= 2 

(Integer  (R/2)) 

R2 

= R1 

+ 

2 

N1 

= BO 

+ 

(B1)(R1) 

+ (B2)(R1)^ 

+ (B3)(R1) 

N2 

= BO 

+ 

(B1)(R2) 

+ (B2)(R2)* 

+ (B3)(R2) 

Using  Linear  Interpolation 
N = N1  + C(N2-Nl)(R3-Rl)/2] 
p = 4.9746N  + (7.9043E-6)N^ 

Phase  IV 

N4  = [6.4375  + .010426T  -( 6 . 8925E-6  )T^  + (4 . 9127E-7)TMM 
F = .1(N4)P 

TAKEOFF  DISTANCE  PROGRAM 

BO  = 13.086  -.00017113A  -( 2 . 0655E-7)A^  + ( 3 . 6861E-11)A’ 
-(2.4156E-15)A'' 
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B1  = -.045635  -(7.8931E-6)A  + ( 3 . 7545E-9 ) A*  -( 9 . 7088E-13)A’ 

+ (6.997E-17)A'* 

B2  = -.001317  -(8.2558E-7)A  + (4 . 0739E-10 ) A*  - ( 8 . 548E-14 ) A* 

+ (5.4964E-18)A'‘ 

B3  = -(1.9097E-5)  + (1.3671E-8)A  -( 9 . 4694E-12 ) A*  + ( 2 . 0434E-15 )A’ 
-(1.4617E-19)A‘‘ 

C = BO  + Bl(B)  + B2(B)^  + B3(B)^ 

If  double  datum  on, 

E = 1.9773  + .56598C 
If  double  datum  off, 

E = .54178  + .65876C 

GO  = -(4.8896E+5)  + (8.4974E+1)G  - ( 5 . 7 35 6E- 3 )G^  + ( 1 . 9 3 7 3E-7 ) G^ 

-(  3.1744E-12)G‘*  + ( 2 . 0446E-17)G® 

G1  = (5.8621E+4)  - ( 1 . 0146E+1 ) G + ( 6 . 8807E-4 ) G^  -( 2 . 292E-8 )G’ 

+ ( 3.7387E-13)G‘‘  -(  2 . 3964E-13  )G® 

H = GO  + GKE) 

If  relative  humidity  < 40%,  K = H 
If  not,  K = 4{[(I-40)/1000]+l} 

LO  = 67.124  + .89509K  + ( 2 . 3306E-5 ) - ( 1 . 6 2 54E-9 ) 

+ (3. 3728E-14)K'’ 

LI  = -9.0995  -(1.0856E-2)K  + ( 2 . 1754E-7 ) - ( 2 . 5 3 2 7E-11 ) K’ 

+ (1.197E-15)K'* 

L2  = (1.4782E-1)  - ( 2 . 166 6E-6 ) K + ( 3 . 4274E-9 ) K*  -( 2 . 7817E-13) K’ 

+ (9 . 3077E-18)K‘’ 

M = LO  + LKL)  + L2(L)^ 

If  winds  calm,  M = K 
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XO  = (4.5704E+1)  + .93429M  + ( 2 . 2265E-5)M*  -( 2 . 338E-9 )M^ 

+ (7.941E-14)M” 

XI  = 7.9472  + .014914M  + ( 9 . 0 70 8E-6 )M*  - ( 7 . 12 35E-10 )M’ 

+ (3.0684E-14)M'* 

X2  = 5.3616  -.0085136M  + ( 3 . 59 14E-6 ) - ( 4 . 59 32E-10 )M’ 

+ (1.9889E-14)M'’ 

X = XO  + Xl(N)  + X2(N)^ 

QO  = 2604.2  -2.1694X  + .0010915X^  - ( 1 . 1119E-':’ ) + ( 3 . 6 6 2E-12  ) X"* 

Q1  = -175.73  + .22601X  -( 7 . 5225E-5)X^  + ( 7 . 7018E-9 )X^ 
-(2.5437E-13)X‘‘ 

Q2  = 2.3549  -.0040102X  + ( 1 . 2 8 32E-6 ) X^  - ( 1 . 3 2 34E-10 ) X^ 

+ (4. 3908E-15)X'* 

I 

Q = QO  + Ql(P)  + Q2(P)* 

50  = -400.79  + 1.5801Q  - ( 2 . 0 2 54E-4 ) Q*  + ( 2 . 4111E-8 )Q^ 

-(8.6737E-13)Q'* 

51  = 16.196  -.024333Q  + ( 9 . 34 84E-6 ) Q ^ - ( 1 . 2 594E-9 ) Q ’ 

+ (4.7522E-14)Q'* 

52  = -.14758  + (2.359E-4)Q  -(1.037E-7)Q^  + ( 1 . 6 016E-11) Q ^ 

-(6.3195E-16)Q'* 

S = SO  + Sl(R)  + S2(R)^ 

MAXIMUM  RANGE  CRUISE  TIME  AND  SPEED 
AT  CONSTANT  ALTITUDE  PROGRAM 

BO  = -1  + (5.0794E-3)H  -(1 . 3968E-3)H*  + (8.254E-5)H’ 
-(1.2698E-6)H‘’ 
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B1  = .05  + .0015159H  + (1.123E-4)H^  - ( 3 . 49 21E-6 )H’ 

+ ( 7.9365E-8)H'* 

N = BO  + BKG) 

BO  = .47803  + .0013417D  + ( 6 . 2 28 7E-6 ) - ( 1 . 626 lE-8 ) D* 

+ (1.6438E-11)D'* 

B1  = .08217  + (4.1209E-4)D  - ( 4 . 5 57 7E-6 ) D*  + ( 1 . 6 777E- 8 ) D ’ 
-(2.001E-11)D‘* 

B2  = (4.2143E-4)  - ( 9 . 4 39 7E- 5 ) D + ( 1 . 2646E-6 ) - ( 4 . 8 5 3 7E-9 ) 

+ ( 5.7222E-12)D‘* 

B3  = -(6.6767E-4)  + ( 8.4671E-6)D  -C  0501E-7)D^  + ( 3 . 6 3 8 2E- 10  ) D ’ 
-(  3.7828E-13)D‘‘ 

M = BO  + BKN)  + B2(N)^  + B3(N)’ 

Ml  = M -C(60-T) (2) (M)/12Q0 ] 

V = (710)(M1-.14)  + 100  -E 

T1  = Dl/V 

FUEL  REQUIRED  FOR  MAXIMUM  RANGE  CRUISE 
AT  CONSTANT  ALTITUDE  PROGRAM 

BO  = 4.54  -.16444A  + .0033932A^  -( 1 . 0283E-4 )A’  + ( 1 . 9 26E-6 ) A” 
-(1.3757E-8)A‘ 

B1  = (3.22E-9)  -(3.6664E-3)A  + ( 8 . 9338E-4 ) A*  - ( 5 . 5 9 39E- 5 ) A’ 

+ ( 1.459  3E-6)A'*  -(  1 . 3281E-8  ) A' 

B2  = (6E-4)  + (1.1203E-4)A  - ( 2 . 3 35 8E-5 ) A*  + ( 1 . 4536E-6 )A’ 
-(3.7144E-8)A'‘  + ( 3.3334E-10)A* 

M = BO  + BKG)  + B2(G)* 
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BO  = -(2. 5399E-3)D  + (9 . 7299E-5)D^  - ( 2 . 3516E-7 ) 

+ (1.4251E-10)D'* 

B1  = 2 + (4.2388E-3)D  + ( 1 . 2 326E- 5 ) - ( 1 . 0 2 9 8E-7 ) D’ 

+ (1.7277E-10)D‘* 

L = BO  + BKN) 

F = L/V 
R = (F)(T)/60 

MAXIMUM  RANGE  CLIMB  AIRSPEED  SCHEDULE 

S = 405.56  -.79075D  + .0011382D^  - ( 4 . 10 18E- 7 ) D ^ 

M = .86  -(2 .1634E-3)D  + ( 7 . 6582E-5)D^  - ( 1 . 1344E-6 ) 

+ (7.2125E-9)D‘*  - ( 2 . 30  3 5E- 11)  D®  + ( 3 . 6 5 8 8E-14  ) D® 
-(2.3062E-17)D^ 

TAKEOFF  AIRSPEED  PROGR.^M 

U1  = 54.023  + (3.4787E-3)G  -( 1 . 9475E-8)G^ 

U = U1  + C(26-P)/2] 

VO  = -1917.1  + 61.604U  -.  70348U^  + .0035661U®  - ( 6 . 6 5 78E-6  ) U** 

VI  = 76.824  -2.4517U  + . 028779U^  - ( 1 . 47  5 3E-4  ) U®  + ( 2 . 7872E-7  )U‘* 
V2  = -.72239  + .023415U  - ( 2 . 79 8E-4 ) U*  + ( 1 . 4 5 9 6E-6 ) U’ 

-(  2.807E-9)U'* 

V3  = VO  + Vl(R)  + V2(R)^ 

MAXIMUM  REFUSAL  SPEED  PROGRAM 

BO  = -43.01  + 6.761G  -.35l59G^  + .0080545G®  -( 6 . 7769E-5  )G'‘ 

B1  = 26.312  -3.8382G  + . 20326G*  -.047022G’  + ( 3 . 9 94E- 5 ) G** 
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B2  = -4.9639  + .72723G  -.038721G*  + (8.985E-4)G’  - ( 7 . 6 38E-6 ) G” 


B3  = .30288  -.044855G  -.0023921G*  - ( 5 . 5 54 9E-5 ) G’ 

+ (4.7217E-7)G“ 

R = BO  + Bl(E)  +B2(E)*  + B3(E)’ 

BO  = -11.412  + 62.185L  -9.0037L^  + .64921L’  -.017455L'* 

B1  = -.2811  -4.2012L  + .703771^  -.058693L*  + .0017461L'* 

M = BO  + BKR) 

OPTIMUM  ENDURANCE  ALTITUDE  PROGRAM 

BO  = 55.333  + .073076D  - ( 9 . 7 8 36E-4 ) + ( 3 . 50 15E-6 ) D ’ 
-(3.9782E-9)D- 

B1  = -1.1  -(8.0597E-3)D  + ( 8 . 0 09 7E- 5 ) - ( 2 . 8 3 36E- 7 ) D’ 

+ ( 3.3032E-10)D‘* 

B2  = (6.6667E-3)  + (1.2541E-4)D  - ( 1 . 40 39E-6 ) 

H = BO  + Bl(G)  + B2(G)^ 

CRUISE  CEILING  PROGR.^M 

BO  = 85.118  -.29117D  + ,00304340^  - ( 1 . 2 8 51E- 5 ) D’  + ( 1 . 6621E-8  )D'* 

B1  = -2.7877  + .025635D  - ( 3 . 306 3E-4 ) + ( 1 . 4 16 2E-6 ) D’ 
-(1.8343E-9)D'* 

B2  = .063327  - ( 8 . 5 289E-4 ) D + ( 1 . 0 8 14E- 5 ) - ( 4 . 6 5 14E-8 ) D’ 

+ (6 .0606E-11)D'* 

B3  = -(6.0468E-4)  + (9.0826E-6)D  - ( 1 . 14 3E-7 )D*  + ( 4 . 9 30UE- 10 ) D’ 
-(6.4567E-13)D'’ 

H = BO  + (Bl)G  + (B2)G*  + (B3)G’ 
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1 PEW  THIS  PPlM-.SHM  ChLCULmTES  THF  FUEL  FLOW  tiWli  LE'KUEL  tlHUTICHL  MILE  FOR  HH 

2 FEM  H-TE  Fl.VlUG  h LOW  I CVEL  MISSION  ONIi  IS  IiEPENDENT  i.iN  4 VHPIhBLES  — 

S PEM  GROSS  WEIGHT. liPMG  COUNT. MOCM  NUNBEP.  HND  TEMPEPHTUPE > C ENT IGPrtDE ' 

10  PRINT  "ENTER  GPO-iS  WT.IiRmG  CT.  MhCN  M.  MNli  TEMFaCENT)  ■ 

11  PRINT 

12  PRINT 

INPUT  O.D.N.T 
•*  ‘ G»G  1000 

N1--0.  3SSl5«0.0O4208l»G 
'?4  r.OSUB  300 
I-O 

•:3  GOSUe  300 

30  S2a3 

70  IF  Si '32  THEN  loO 
•0  S«S2 
GOTO  300 
100  I«l 
1 1 0 Gosue  300 
120  S3-S 

no  IF  SI  - S3  THEN  300 

140  S3-S3 

If.O  l»I*l 

130  GOSUB  300 

l.'O  GOTO  130 

300  Il»(Sl-S3  ■ ■n3-S3> 

310  Ml»M 

330  I»I-l+Il 

331  I»INT(I> 

333  GOSUB  300 

333  S2«S 

334  t«I+l 

33?  GOSUB  300 
333  S3«S 

337  S»S2+' I l»vS3-S3  ' ' 

340  GOTO  300 
3S?  PRINT 
3S3  PRINT 

5U0  R«S*3*k4.  3733E-03+0.037743*D  >-»N»‘3 

301  R3»R 

303  Rl«34lNT(R  3.> 

304  R3'«Rl*3 
303  J«l 

305  IF  J«2  THEN  311 
30'?  P»RI 

3!0  GOTO  31S 

311  R»R2 

3 1 ■?  B0-?  .335  3-1.  SSS-R  ♦ 3 . 0353  *R  f3- 1 . 073  1 »R  r 3*0 . 1 7375*R  r4 
330  e0»B0-0. 01  30S5»Pt?(.3.  '52(5E-04*Ptii 

3 JO  B 1*305.  3013-34S.'?317*P«-St  .3r3355-pr3-15.5531S»Rf:3<-l  .334  4 33*pt4 
340  et»ei-0.03S5  33  3*RT5»3.SS3  3S?E-04+RTi!i: 

■i50  B3*-1053.  133»1331 . 34»R-487. 4333*RT3»'?1 . 3533 “R  ^ 3-S.  333S33*Pt4»0.  5?53S74*P  ^5 
330  B2>*B3-0.003?05?35*PT3 

570  B3-1300.  l43-r-‘50.  1 3S*R+7SS.  85 1 3>Rt3-i53. 5733*R  " 3>  1 5. 038  n*R  r4 
530  B3«B3-0.  73741  3-^+F  r5+0.0l3707»RtG 

3'50  64«-S34. 3375+ 1000. 44 3*P-408.  745 1 »R  T3-I-S0. 0S31 4*R t 3-3.  o:s‘iS*RT4 

400  B4«e4*0. 3SS3537*Pt'5-7.  7303 1 7E-03*-R T3 

430  N»B0^Bl*M+B3*Mf3+B3»MT3»64*MT4 

440  IF  J«3  THEN  480 

450  Nl«N 


4??  Ja2 
4(50  GOTO  Sit 
4^5  Ra2 
4S0  N2aH 

4>0  HatlU' M2-H1 '"'(PS-PP  2 

r.OO  PEN  COMPLETED  C MLCULhT  I nil  OK  IMTEPMEDIOTE  » BY  LltlEOP  1 tITEPPOLOT  ION 

? 1 0 P«4  . '4t>»M«-  C . 'f'04  'E  OB-cN  f2 

'■20  N4a('j.4  5r'5  + U.01042G*T-E.:S025E-0'i-»Tt2  + 4.  ?12rE-Or*Tt-0  “M 
r?50  F*(0.  1 *M4»P>-rlOOO 
F-INT.  P ' 

'iJO  PPINT  "GROS';  NT»n>1000 
■541  PRINT  ■•T'5  = "'5  DC’ 'D''Ma”M 
?42  PRINT  "TEMPa'T 
?4J  PRINT  "PEF  H’-RO 
?44  PRINT  "Na"N 
^4-5  PPINT  '•LeFUEI.  NMa"P 
ri^O  PPINT  "FLiEL  FLOW*  "F 
^51  PRINT 
nr?  GOTO  10 

GOO  B0“22.  or?-  31  . r?4  > I *-4  I . ?3»  P2-n.  ons^r  I f 5 
G 10  Bl  *-  154.  -^0  + 21  r.  5 1 ' I -2t.  1 . 23*  I f 2 co? . '^OS* I TS 
G20  B2’405. 00-525.  'G*I  *•002.4'^-  I f2-0:3.  r37«I  T3 
6 30  B3*-445. 62»542.  OO'  I -6  I 1 . 55»  n2  + '^2.  :5'?4*  1 1 3 
^40  B4*  104. 70-204. 42  ♦ I ••225. 0'3->  I '“2-  35.  lOO-*  I TS 
K50  S*B0+61-»M1  (•62^01  »2*-B3»ni  t-3+e4»Ml  M 
660  RETURN 

COO  00 a -2.  3207-0. 26  3 1 6 »D‘0. 007  3.527’D '2-'.  5 1 3E-05»DT  3 + 3. 5 306E-07*D  »‘4 
010  00*00-7.  70E-lO*Di5'-6. 4624E-  1 3-Df6 

320  HI  *4. 0 35  + 1 . 0'^5'>*D-o.  0 30657  -Df  2 + 3.  1 012E-O4*Dt  3-1 . 5276E-06  + DT4 
3 30  0 1 *0 1 + 3 . 400E -00 - D 1 5- 2 . 06 026 - 1 2 ♦ D 1 6 

40  (02*10.  204-1 . 071  '^^D+O.  0310'74-D*2-3. 2073E-04*D»'  3 + 1 . 5'?5E-06’Dt4 
050  »02*h2-  3. 60OOE-OO  ••DT5  + 3. 06  34e-l2*DT6 
060  0 1 *mO  + h 1 'M I +02 *H  1 12 
070  RETURN 
000  END 
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List  of  Variables  for  Program  1 


Drag  count 
Temperature  (°C) 


Mach  number 


Result  of  lower  graph,  Figure  B1 


Guidelines,  numbered  top  to  bottom  consecutively 


Transfer  Scale  calculated  as  function  of  I 


Transfer  Scale  calculated  as  function  of  D 


Transfer  Scale  calculated  for  upper  guideline 
Transfer  Scale  calculated  for  lower  guideline 


Relative  Transfer  Scale  location  between  guidelines 


Reference  number 


Even  reference  number  below  actual  reference  number 


Even  reference  number  above  actual  reference  number 


Result  of  lower  graph,  Figure  B3 
Result  of  lower  graph,  Figure  B3  for  R1 
Result  of  lower  graph.  Figure  B3  for  R2 


Result  of  lower  graph.  Figure  B4 


Pounds  of  fuel  per  nautical  mile 


Fuel  flow 


Program  2 


I PEM  THIS  PPi.ir.RtiM  CHLCHLftTES  THE  Tfil  EuFF  DISTHHUE  REOUIPCO  FOR  HH  fl-7E 
PEM  IT  IS  PEPEtlDEHT  ON  ■?  VRKUTBLES  — 

J.  PEM  GROSS  HEIGHT, PNir,  HLT  I riJPE,  TEMP.  DRfiG  CO'JUr.PfLHTIVE  HI.IMII1I  TV.  WINPS 
•i  REM  RNWY  SLOPE, I. ENTER  OF  GPHVITV  LOCHT ION.  FLAPS.  ANH  DOUBLE  DATUM  STATUS 
PRINT  "INPUT  ALT.  TEMP.  DC. GH" 

10  INPUT  A.e.D.G 

II  1=50 
12  L=10 
IS  N“t 
14  P-27 
20  R-25 

too  eO=l3.0SS-O.000irn3.^A-2.OSb5£-O.’.At2+3.6'StiIE-I  I*Af3 
im  eO-eO-2.4l5E.E-15-At4 

1 10  ei— 0.045S35-7.S031E-0S.H. 3. 7545E-0S*At2 

1 1 1 Sl-ei-S.  70SSE-t  3*A  t 3*S.S07E-l7»A*'4 

120  e2--0.00t3l7-S.255SE-0r+A+4.07S'?E-10*At2 

121  e2=B2-S.54SE-14'Mt3.5.4'3S4E-lS»At'4 

I 50  B3=-l  . SOSrE-OS+l . 3SriE-034A-S.4S'?4E-12*At'2+2.04  34E-l5-*At;3 

140  B3-B3-l.4otrE-lS4Flt4 

ISO  i:=B0+Bl*B.'B2+B42^B3»'Bt3 

ISO  IF  D-1  THEN  ISO 

170  E=0. 54 1 rS+O . S5S7S4C 

ISO  GOTO  200 

1=0  E-1.S773+0.5S5SS4C 

200  G0»-4.SSSSEh5Sh3.4S74E+01.'G-5.  7S5SE-034G f'24 1 . S'37'3E-07+G  f 3- 3.  1744E-12»Gr4 
-to  G0-G0+2.044SE-17'■G^5 

220  Gl=5,SS21E+04-1.0l4SE+0l  '.G.-S.SS07E-04*Gl'2-  = .2S2E-03*Gf  3i-3.7387E-l3*G»'4 
230  Gl»Gl-2.3'?S4E-l!34G55 
-40  H-G0+G1»E 
250  j-0 

2S0  IF  I<40  THEN  2S0 
270  J«<I-40>  1000 
2S0  K-H4J+H 
235  IF  L-0  THEN  340 

S-^O  LO-S.  7124E>01  rS.  S50SE-01  *-2.  330SE-05*Kt2- 1 . S254E-0S4K  r3+3.  372SE-14*K  f4 

500  L1--S.09S5-1 . 0SSSE-02-K  ^2.  1 754E-074K >2-2. 5327E- 1 1 -.P  53+ 1 . r?7E-15*K54 
: 10  L2-1.4732E-01-2.  1S>3SE-0S+P +3. 4274E-0S4K52-2. 73 1 7E- 1 3+K  5 3+S. 3077E- 13*Kt4 
320  M-L0+H+L+L2+L52 
-.30  GOTO  350 
540  M-K 

350  ;;0-4.5704E+01+'?.  342'?E-0l+M*2. 22S5E-05+Nt-2-2.  333E-0-5+N5  3 + 7,  S4  IE- 14+M54 
3F.0  Xl-7.  S472  + 1 . 4'5  14E-02*N  + '4.0703£-0S'M52-7.  1 235E- 1 O+M  + 3 + 3.  0SS4E-14+N54 
370  ::2-5. 5elS-S.513SE-03*N+3.5?14E-0S+H52-4.5S32E-10+N53+l .S3SSE-14+Mt4 
.'SO  .:-:<0+Xl +N+X2+NT2 

3 50  00-2.  S042E+03-2.  lS'?4+;<+ 1 . 0'5l5E-03+Xt‘2- 1 . 1 1 r5E-07*X+3+3.  SS2E- 12+X M 
400  0I--1 . 7573E+02+2. 2S01E-01  +X-7. 5225E-05+:;t'2+7. 7013E-0'?*.;t  3-2,5437E-13+XT4 
410  02-2.354'5-4.0102E-03+X+1.2332E-OS+X52-l.3234E- lO+Xi J*4. 3S03£-15+Xt4 
420  0-00+0 1 ♦P+02+PT2 

■t  ^o  S0--4. 007'=£+02  + l . 5.'?0I  +0-2.  0254E-04+0T2+2.  41 1 1 E -00 *0  + 3-3.  S737E-13+0+4 
440  Sl-l.SlSSE+01-2.4  33  3E-02*0+S.3434E-0S+Q+2-1.25  54E-0'5+0t3+4.7522E-14*C!t4 
450  S2— 1.4753E-01  +2.  353E-04+0- 1 . 037E-07+U +2+ 1 . SO  1 SE- 1 1 +0t  3-S.  3l‘55E- lS»Qt4 
4S0  S-S0+Sl+R+S2+Pf2 
470  S-IMTkS) 

475  PRINT  " FOR" 

430  PRINT  "GW-"G"  ALT»"A"  TEMP»"B"  DC«"D"PH»" I"HDWD-"L 

432  PRINT  "RNWV  SLP»"N"’.  CEN  GRAV-"P"FLAPS*"R 

433  PRINT 

530  PRINT  "TAKEOFF  POLL  DIST«"S 

531  GOTO  3 

532  END 
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List  of  Variables  for  Program  2 


Variable  Definition 


Runway  Altitude  (feet) 


Double  datum  status  (1  indicates  "with") 


Gross  weight  (lbs.) 
Relative  humidity  (%) 


Headwind  (kts.) 


Runway  slope  (%) 


Center  of  gravity  (%) 

Flap  position  (degrees) 

Result  of  upper  graph,  Figure  B5 


Unadjusted  ground  roll  distance.  Figure  B6 


Adjustment  factor  due  to  relative  humidity 

Ground  roll  distance  (GRD)  adjusted  for  relative 
humidity 


GRD  adjusted  for  runway  slope 

GRD  adjusted  for  the  center  of  gravity  location 
True  GRD  (also  adjusted  for  flap  position) 


1 REM  THI5  FRor.RHii  the  li-rE  pmnge’ hif-ipeed  rud 

2 REM  TIME  OF  FLIGHT  mHD  VI  DEPEMDEHT  Oil  o ''RRIHBLE?  -- 

3 REM  GROSS  HEIGHT . FiLT  I TUDE • DPhG  COUN  1 . TEHPEPRTIJPE • U I MDS • HMD  DIPTHNCE 
S PRItir  "INPUT  GW- RLT.  DC  - TEMP- tC  ■ !HDUD«  tUSTFiNCE" 

10  INPUT  G.H.D'T.L-DI 
30  G^G.  1000 
40  HsH  1000 

50  FiO=-l  +5. 0r34E-O3»H-  ! . SOESE-OS+H  ^3-3. 4Sl3E-0E*H^3<-r.  0 3E5E  -0S»H^4 
EO  Fll=0. 05+0. 0015153'^H'- 1 . 1 33E-04-»H  ^3-3 . 403 1 E-OE-^Ht  S + T.  '?3E5£-03»H’'4 
rO  N=R0  + R1-»G 

SO  e0=0. 4 7303-0. 00  1 34 1 r»D  + E.  33SrE-0E-D  > 3 - 1 . E3E 1 E-OS^^D  ^3f  1 . E4 SSE- 1 1 - D <'4 
35  81=0.03317+4. 1303E-04  + D-4. 5577E-0E  + Dt2+1 . E777E-03  + D t'3-i , 00 1 E- 1 1 *D f 4 
00  B3=4.3143E-04-0.4  307E-05-D+l . 3E4EE-0E  + D t'3-4 . ;3537E-00-D f 3+5 . 7333E-  13  + Df4 
05  B3  = -E. E7E7E-04  + S'. 4E7 1 E-OE  + D- 1 . 050 1 E-07  + DT3  + 3. E333E- 1 0 + D rS- 3 . 7333E- 1 3*D 1 4 
100  N=60+B1 +N+B3*Nr3+B3+Nr3 
110  M = M-',  ' >EO-T  •+3-'N  • 10+130>  ' 

130  V=710+' M-0. 14 - + 100-L 
130  T1=D1  V 
135  V=INTkV> 

140  PRINT 
150  PRINT  "FOP" 

lEO  PRINT  ■•Gl<i=-G"  FtLT  = "H"  00  = " D" TENP  = " T " HIiND  = "L" D IST  = " D 1 

lEl  PRINT 

1E3  PRINT  "GROUND  SPEED="V"  TINE  OF  FLIGHT="T1 
170  END 
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Variable 

G 

H 

D 

T 


D1 

N 

M 

V 

T1 
AO,  BO 


A1,B1 


B3 


List  of  Variables  for  Program  3 

Definition 

Gross  weight  (lbs.) 

Altitude  (ft.) 

Drag  count 
Temperature  (°C) 

Headwind  (kts.) 

Distance  to  fly 

Result  of  first  chart,  Figure  B8 

Cruise  Mach  number  (adjusted  and  unadjusted  for  T) 
Ground  speed  (kts.) 

Time  of  flight 
Coefficients 
Coefficients 
Coefficient 
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Program  4 


10  PEM  THI?  PPOr.PflM  CRLCULPlTEi  FUEl.  PEOUIPEI'  POP  Mti:'  PrtNGE  OT  CONSTPINT 

11  PEM  mLTITUDE  pop  htl  Pl-TE  hMP  I j DEPEUDENT  OH  p VmPIOBLES  -- 

12  PEM  GROSS  WEIGHT. OLriTUDE.DFOG  OOUMT-TFUE  OIPSPEED-  mHD  T IME ' M IKUTES > 

20  PRINT  "ENTER  nPOSS  WT.iiLT.DRmG  CT  . ThS  • T IME  M I NOTES " 

2 I PRINT 
22  PRINT 

?0  INPUT  G.M.D.V.T 
2^  PRINT  "GROSS  WT=i"G 
print  "mltitude= 'm 
57  PRINT  "DRMG  i.Ol.lNr.s  D 
5S  PRINT  “TRUE  hIRSPEED»"V 
?•..  PRINT  TIME  OF  PLIGHT*“T 
hO  G’G  1000 
‘.‘i  M = M lO'iO 

0 0 eO»4.  '4-0.  1P44.*  'H+o.  0055022  >0^2  ■ 1 . 02S 5E-04 ♦pt-j 1 1 . ■52SE-0S i-M 1 4-  i . Sr^PE-OSiMP? 
'<■'  et^-S.  -.'SPJf:  ..t  Vhi-S.  ?55SE-t>4»R^2-?.  5'?5'^E-05*M'T5*- ! . 45?.5£-0<5 i'rtt4 - ! . 52S 1 E-t3.S ♦•M  T*! 
.0  B2*r;E-04*i . 12O5E-04  >h-2.  555SE-0S->M T2+ 1 . 455SE-0S-M T5-5 . 7 1 44E-0Si-M t'4 
E:2*e2»  5.  5554E-10*Mf^i 
'*0  N = B0+6 1 *G*b2  *G  * 2 

100  M03-2.  ? 5';-4E-0.  »D»'4.  72'5^E-05-D  t2-2. 551bE-074DT  2*l . 4251E-104DT4 
1 10  Ml=2*4.25SSfe  0:iii+ 1 . 222bE-0S-Df2-l . 029SE-07  + D f 3+ 1 . 7277E-104DT4 
120  L=»M0+M1*N 
150  PaL'V 

140  p»F4T  so 

141  La^  INT. L»  1000  ' ■ IOhO 

142  PaINT' F^ 

145  RaINT' R ' 

147  PRINT 
14:5  PRINT 

150  PRINT  "LBFUEL  NMa"L"FUEL  PLOWa'F 

1?5  print  "FUEL  REOUIRED="R 

15S  PRINT 

157  PRINT 

ISO  GOTO  20 

170  END 
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D Drag  count 

V True  airspeed  (kts.) 


Program  5 


1 PEM  THIS  PPOGRPIH  C flLi. IJLI=ITE3  THE  CLIMB  filPCPEED  OF  HM  m-TE 
;■  PEN  ' IHrUCMTEIi  HIRBPEED  BELOW  20.000’  ' 
j REM  ' MHCH  HIJH6ER  mBOVE  2O.OO0’> 

1 0 D=*0 

12  PPIHT  "CLIMB  HIPBPEED  SCHEDULE" 

15  PRINT  "DRFii:  CT  CLIMB  AIRSPEED  CLIMB  MACH" 

IE  PRINT  '■  ( IhS  to  20000’''  'ABOVE  20000’'" 

20  S =403. 56-0.  f'S075-''D+0.  Out  1332*DT2-4.  10 1 3E-07'fD T3 

21  S=INT(S.' 

:0  M=0.:56-2.  lE34E-03''D.-7. 6502E-05-^Dt2-l . 1 344E-0E'-DT3  + 7. 2 1 2:'E-0S*D M-2 . 3035E-U 
40  M=M  + 3.o533E-l4  4DtE-2.3062E-17*Dt'7 
42  M=M»100e 
44  MsINTCMI' 

46  M=M.  1000 
55  PRINT  D-S.M 
CO  D=D+30 

^■O  IF  D<310  THEN  20 
30  END 
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D Drag  count 

M Mach  number 

S Calibrated  airspeed  (kts.) 


31 


Program  6 


400  PEN  THI'i  PROCPhM  CPiLruLRTES  THE  ThKEOFF  hIR'JPEED  OF  hri  h-7E 

401  PEN  UNDER  VRPVING  OROSO  UEIOHTS-  FLliP  POS I T luNi*' 

402  PEN  hHD  L enter  OF  GRhV I T','  LOCRT ION'S 
4'?;5  R=20 

4'?'?  P = 20 

500  0=20000 

501  PRINT  "FOP  GROSS  WEIGHT="G 

502  PRINT 

503  PRINT 

504  PRINT  "FLP.PS  CG  TAKEOFF  AIRSPEED" 

530  Ul=5. 4023E+01+3. 47S7E-03-^G- 1 . •?475E-0S*G  f 2 

540  U=U1  + ' 2t'-P.’  2 

550  V0=:- 1 . 3 1 7 lE  + 0 i-rP. . 1 P04E  +0 1 'rU-7 . 034SE-0 1 -Ut2  + 3 . 5P6  lE-OS'^U  t 3-P . P57SE-0P*U  f 4 
530  VI  =7.  PS24E  + 01-2. 4517-U  + 2.  S77'?E-02-Uf-2-l . 4 753E-04^U  f 3 + 2.  7S72E-07+Ut4 
570  V2=-7.  2233E -01  +2. 341 5E -02  + 1.1-2.  7'3SE-04  + UT-2+  1 . 45:3PE-0b+U ^ 3-2 . S07E-03+U 1 4 
5 y 0 V S = V U + ' .■'  1 - R + ’’.•'  4 + R’  + 2 
590  V4=INT' V:::  .' 

600  PRINT  R,F,V4 
610  R=R+5 

620  IF  R 40  THEN  PiO 
625  GOTO  530 
6301  P = P + -J 
631  R=20 

635  IF  P 35  THEN  P50 
640  GOTO  530 
P50  G=G+300Ci 

651  R.=20 

652  P=20 

655  PRINT 

656  PRINT 

657  PRINT 

6601  F’RINT  "FuP  GPOS'5  W£IGHT  = "G 

662  PRINT 

663  PRINT 

669  IF  G 42000  I HEN  710 

670  GOTO  530 
710  END 
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List  of  Variables  for  Program  6 


Variable  Definition 


Flap  position  (degrees) 

Center  of  gravity  (%) 

Gross  weight  (lbs.) 

Unadjusted  takeoff  airspeed 

Takeoff  airspeed  adjusted  for  center  of  gravity 
Actual  takeoff  airspeed  (adjusted  for  flap  position) 


VO, VI, 
V2  ,V3 


Coefficients 


Program  7 


1 KEN  THI?  PKOi'.RnN  CflLCULflTES.  THE  MNNIMUN  REFUSAL  SPEED 
L REN  FOR  AN  m-;‘E  iI-IHi:  ANT  I -skid 
S REN  IT  I'o.  DEPENDENT  ON  S VARIABLES  -- 

4 REN  GRO-jS  UEIGHT.  TENP.RNN')'  LENGTH-KNNV  ALTITUDE*  AND  DOUBLE  DATUM  STATUS 

4 PRINT  -INPUT  ALT. TEMP. RNWV  LTH. GW. DOUBLE  DATUM" 

10  INPUT  A.B.L.G.D 

IS  G«G  1000 
SO  L=«L  1000 
~0  PRINT  ■•ALTITUDE=»”A 
PRINT  "TENP-'B 
"S  PRINT  "Rrun'  LTH»"L*1000 
~5  PRINT  "GROSS  WT«"G»1000 
'T  PRINT  "DD»"P 

100  BO=l >. OSE-0. 0001 '1 1 ?-A-S.OGSSE-Or>ATS+S. ESBlE-1 1 »A^3 

101  BO-BO-S. 4lSBE-tS»At4 

1 10  B1  «-0. 04SBSS-T.  AA31E-06*A*  ?.  .'S4SE-0'4*A T2 
III  ei«ei -4. 'ossE-i 3*At 5»B.4orE-i r»At4 

ilo  es*-o. 001 31  .’-o. ssssE-or*A*4. o^?•:>E-^o♦A^s 

l-l  5S*eS-;?.S4SE-l4*At  ? t-S . 4SG4E- 1 S*A  f 4 

ISO  e?«-l . POorE-OS*-!  . SSriE-OStA-'?.  4i;44E-  IS^AtS  + S.  04S4E-15*At3 
140  B3-63-1 . 4G1  :E-14-.Af4 
ISO  C^BO+ei+e+BS'BfS+BS^Bf? 

IcO  IF  D«1  THEN  ISO 
i"o  E’O. S4 1 r'.?*o. sssrA»i' 

ISO  GOTO  SOO 

140  £»1 . 4,'73«.o.  SSS4;3,i; 

.;0o  B0»-4S.  01  ♦£.  TGI  .G-0.  SS  IS4>|- m).  00SO54S*Gt' '-G.  ~rG4E-0S*G’^4 

5 to  Bl  *2G.  312  - S3S2  ‘Gt-O.  20  32G  *0 1 2-0 . 004 "022*G  f 3 + 5. 4 '>4E-os.G  t-4 
220  B2«-4.  SGSO+O.  r2r22»G-0. 0 30221  H;t2»S.  4S5E-04-GL5-r.  g:'SE-0G»GK4 

- ■'  B3’0.  302S0-0. 0440SS«G»0. 0023321  ^G r2-S . SS44E-0S-.G ^ 4*4 . 721  rE-or»G*'4 
S-u;.  PaB0  + Bl*£*B2.-Ef2*6.3*Et'3 

^'O  B0»  - I 1 . 4 12»G2.  10S*L-'?.  0037*L  I'j^O.  G4'321  *L  f 3-0. 01  74S5»L  '■4 
^GO  ei  «-0. 2S1  I -4. 20i;*L<-0.  70  3r7*Li2-0.  OS0G'3  3*Lf3  + O.  OOI  '4G1  *L  f4 
370  M«60+ei«R 

_.'l  N»INT'M' 

2 "2  PRINT 

2 "s  PRINT  " Takeoff  factor  »“E 

2 V print 

7 m print  "MA:,  refusal  speed  « '.M 

201  PRINT 
^40  GOTO  4 
.00  END 


List  of  Variables  for  Program  7 


Variable 

A 

B 

L 

G 

D 

C 

E 

R 

M 

BO ,B1, 
B2,B3 


Definition 

Runway  Altitude  (ft.) 

Temperature  (®C) 

Runway  length  (ft.) 

Gross  weight  (lbs.) 

Double  datum  status  (1  indicates  "with") 
Result  of  upper  chart.  Figure  B5 
Takeoff  factor 

Result  of  first  chart,  Figure  B12 
Maximum  refusal  speed  (kts.) 

Coefficients 


85 


F 


Program  8 


l PEM  THIS  PPOCRHM  CHLCUl HTES  THE  OPTIMUM  EMDURHNCE  RLTITUDE 
^ REM  OF  AN  rt-TE  HT  VHP'i'IHG  GROSS  WEIGHTS  HMD  DPHG  COUNTS 

4 DIM  BC3] 

5 G»IS 


10  PRINT  •OPTIMUM  EMDUPHMCE  ALT  " 

30  PRINT  "GROSS  WT  DRHG  CT  OPT  END  HLT'' 

50 

so  ec  3 >5?. 333+0. 0r30?S»D-S. 7S3SE-04  + D+3+3. SOl^E-OS+Df "-S. 57S3E-0‘aTDt4 
'.■»0  Bt  1 >-l.  l-S.O5'?7E-O3+D<S.O0SrE-O5»-Dr2-3.SS3SE-67  + Df3  + 3.  i032E-io*D+4 
100  et  3 ]»S.SS6  7E-0S  + l . 354  1 E-04  + D- 1 . 403'?E-0S*Dt2+5. 3033E-0'=>‘DtS-t;.  6'^  ISE- 1 ^ + I 
110  H=BC3]+BC1]»-G+BC3}*GT3 
1 15  r»INT. H+1000> 

113  X=G*1000 
ir?  PRINT  X.DjZ 

130  D=D+30 

131  IF  DOlO  THEN  80 
123  D'O 

133  IF  G<45  THEN  50 
140  END 


Program  9 


1 REM  THIS  PR0^■•.R^TM  CHLCULH TES  THE  CRUISE  CEILING  OF  HN  H-TE 

2 REM  UNDER  VHRVING  CROSS  WEIGHTS  NND  DRRG  COUNTS 

4 DIM  BC4] 

5 G=*19 
ij  D“0 

10  PRINT  "CRUISE  CEILING" 

20  PRINT  "GROSS  WT  DRAG  CT  CRUISE  CEILING" 

50  G^G+3 

SO  BC4]=S5.  1 10-0.2'=' U 7 ‘■D^'O.  0030434  ♦DfS- 1.235  lE-05-*DrB+l.SB21E-08»Dr4 
90  BC 1 ]=-2.  7377  H>. 0253 35 »D-3. 303 3E-04»D 12+ 1 . 4 132E-03  + D 13- 1 . 334SE-0'a+Di4 
100  BC  2 ]“0. 033327-3.  5239E -04  <-D+ 1 . 03 1 4E-05  + D12-4 . 351 4E-03  + D 13+3. 0303E- 1 1 *Dt4 
105  BC  3 I=-3. 0433E-04*'?.0323E-03»D-1 . 14  3E'07+Di2+4.  '3304E- lO  + D i3-3. 4537E- 1 3*D14 
I 10  H-BC4  1+BC  1 I+G  + BC2I+G12  + BC3I+G13 
1 15  2*INT<'H+1000> 

113  X=G*1000 
11?  PRINT  X»Dj2 

120  D»D+:30 

121  IF  D<310  THEN  30 

122  D=0 

123  IF  G-45  THEN  50 
140  END 


Variable 

G 

D 

H 

Z 

X 

B3,B4 


List  of  Variables  for  Program  9 
Definition 

Gross  weight  (lbs.  times  lOOC) 

Drag  count 

Cruise  ceiling  (ft.) 

Cruise  ceiling  (integer  format) 

Gross  weight  (lbs.) 

Coefficients 


k. 


89 
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APPENDIX  E 


USER  INFORMATION  FOR  PROGRAM  1 


Number  of  Steps:  1336 


Computation  Time:  90-110  seconds 


ENTER 


PRESS  KEY  DISPLAY 


gross  weight/1000 


drag  count 


drag  count 


mach  number 


mach  number 


Unusable  number 


drag  count 


mach  number 


mach  number 


lb . fuel/nautical  mile 


fuel  flow 
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0 3 

2* 
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07 
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04 
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4 5 

43 

65 
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STO 

43 
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09 
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45 
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4 3 
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05 

5 

366 

06 

06 

35 

■t- 

367 

> ■> 

IHV 

03 

5 

363 

G£ 

93 

369 

04 

04 

00 

0 

370 

01 

01 

06 

6 

1 

V < i 

43 

RCL 

03 

■% 

1>T  2. 

09 

0'5 

04 

4 

*■  '! 

• > 

stn 

52 

££ 

-5  r 4 

0:5 

03 

94 

?r*5 

* ^ 

\ 

RCL 

01 

1 

5 

07 

07 

02 

*■ 

65 

:< 

373 

oi 

1 

4 3 
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